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We consider the antiferromagnetic Heisenberg and the repulsive Hubbard model on two iV-site 
one-dimensional lattices, which support dispersionless one-particle states corresponding to localized 
states on triangular trapping cells. We calculate the degeneracy of the ground states in the sub- 
spaces with n < n max , n max oc N magnons or electrons as well as the contribution of these states 
(independent localized states) to thermodynamic quantities. Moreover, we discuss another class 
of low-lying eigenstates (so-called interacting localized states) and calculate their contribution to 
the partition function. We also discuss the effect of extra interactions, which lift the degeneracy 
present due to the chirality of the localized states on triangles. The localized states set an extra 
low-energy scale in the system and lead to a nonzero residual ground-state entropy and to one (or 
more) additional low-temperature peak(s) in the specific heat. Low-energy degrees of freedom in 
the presence of perturbations removing degeneracy owing to the chirality can be described in terms 
of an effective (pseudo)spin-l/2 transverse XX chain. 
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I. INTRODUCTION 

The thermodynamics of strongly correlated lattice 
models is generally unknown. Although analytical 
(like conventional Green's function technique, dynamical 
mean-held theory etc) and numerical (like series expan- 
sions, quantum Monte Carlo algorithms, density-matrix 
rcnormalization group algorithms etc) methods being ap- 
plied appropriately in particular cases may yield desired 
thermodynamic characteristics with required accuracy, 
seeking for new approaches permanently attracts much 
attention of theoreticians. One interesting idea for calcu- 
lating thermodynamic quantities for strongly correlated 
systems, which is related to the concept of localized one- 
particle states^ has been suggested recently for some 
spir>2^ and electron^ models. The localized nature of 
one-particle states for certain classes of lattices allows to 
construct exactly the relevant many-particle states and 
to estimate their contribution to thermodynamics using 
classical lattice-gas models which are much easier to in- 
vestigate than the initial quantum many-body models. 

In previous investigations of localized states performed 
for the antiferromagnetic Heisenberg and the repulsive 
Hubbard model on highly frustrated lattices^— mainly 
bipartite trapping cells (e.g., a single bond or equilat- 
eral even polygons) were considered. These cells have a 
nondegenerate ground state for the one-particle problem. 
Here we extend the discussion of localized states to highly 
frustrated lattices with non-bipartite triangular trapping 
cells. A non-bipartite cell may have a degenerate ground 
state for the one-particle problem that may lead to new 
effects. For example, an equilateral triangle has a two- 
fold degenerate one-particle ground state, which can be 



related to the chirality degrees of freedom associated to 
a triangle, see, e.g., Ref. [g. To be specific, we consider 
(i) a one-dimensional (ID) lattice which consists of cor- 
ner sharing "double-tetrahedra" (the double-tetrahedra 
chain) and (ii) a frustrated (cylindrical) three-leg ladder 
having a triangular arrangement of rungs (the frustrated 
triangular tube), cf. Fig.[TJ Both lattice geometries have 
been considered in the literature, see, e.g., Refs. for 
the double-tetrahedra chain and Refs. [11H18I for the tri- 
angular tube. Note that triangular-tube geometry is re- 
alized for the copper spins in [(CuCl2tachH)3Cl]Cl2 i 14116 

In what follows we consider two concrete models of 
strongly correlated systems on these lattices, namely 
the spin-1/2 Heisenberg model and the Hubbard model, 
and discuss the consequences of the localizcd-magnon or 
localized-clcctron states in combination with the addi- 
tional chirality degrees of freedom. We mention that 
some similar ideas have been elaborated recently for the 
Hubbard model on decorated lattices^ for the coupled 
tctrahedral Heisenberg chain 7 -^ as well as for the frus- 
trated Heisenberg spin tubeJ^ However, in these refer- 
ences the concept of localized states has not been used 
to discuss low-temperature thermodynamics for thermo- 
dynamically large systems. Furthermore note that some 
preliminary results of our study were announced in a con- 
ference paper J£ 

The paper is organized as follows. In Sec. |TT]we discuss 
the one-particle spectra of the spin and electron models. 
In Sec. IIIII we briefly illustrate the construction of in- 
dependent localized many-particle states and calculate 
the contribution of these states to thermodynamic quan- 
tities. Then, in Sec. IIV1 we illustrate how we can go 
beyond the independent localized states taking into ac- 
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FIG. 1: (Color online) Two ID frustrated lattices with tri- 
angular trapping cells: (a) the double-tetrahedra chain and 
(b) the frustrated (cylindrical) three-leg ladder (or frustrated 
triangular tube). The exchange or hopping integrals along 
the equilateral triangles are J2 > or ti > (bold bonds) 
whereas all other exchange or hopping integrals are Ji > or 
ii > (thin bonds). 



count additional low-energy excitations. Finally, in Sec. 
|V]we consider symmetry-breaking interactions which lift 
the degeneracy related to the chirality. We end up with 
a summary of our findings in Sec. IVII 



II. HEISENBERG AND HUBBARD MODELS 
ON ONE-DIMENSIONAL LATTICES WITH 
TRIANGULAR TRAPS 

In our study we consider the two ID lattices shown in 
Fig. [T] The double-tetrahedra chain [panel (a)] may be 
viewed as a generalization of the diamond chain (the ver- 
tical bond in the diamond chain is replaced by the equi- 
lateral triangle). The frustrated triangular tube [panel 
(b)] may be viewed as a generalization of the frustrated 
two-leg ladder (again the vertical bond is replaced by 
the equilateral triangle). The essential geometrical ele- 
ment of the considered lattices are these equilateral tri- 
angles (which act as trapping cells, see below) together 
with the surrounding (connecting) bonds attached to the 
sites of these equilateral triangles. In order that the con- 
necting bonds should prevent the escape of the localized 
magnon (electron) from the triangular trap, each bond 
of the trapping cell together with two of the connect- 
ing bonds attached to this trapping-cell bond must form 
an isosceles triangle, i.e., the two connecting bonds must 
be equal to each other. As a result, the considered lat- 
tices, owing to destructive quantum interference, support 
localized one-particle states. We note here that the lat- 
tices with triangular trapping cells may be constructed 
in higher dimensions too, see Refs.[^ and[2p|. 

On these ID lattices we consider the spin-1/2 Heisen- 
berg antiferromagnet with the Hamiltonian 

IL ^./,,s, -s, hS z , S Z =Y. S * 

(ij) 



and the repulsive Hubbard model 

H c = H o,<r + U Y1 rc ijt n u , U > 0, 
<r=t4 » 

(ij) ' 

We use standard notations in Eqs. (|2.ip and (|2.2|) and 
imply periodic boundary conditions. The exchange or 
hopping integrals acquire two values: J 2 > or ti > 
along the equilateral triangles (bold bonds in Fig. [T]) and 
Ji > or t\ > along all other bonds (thin bonds in 
Fig. [1]). It is convenient to label the lattice sites by a 
pair of indeces, where the first number enumerates the 
cells (m = 1, . . . , Af, Af = N /4 for the double-tetrahedra 
chain, Af = N/3 for the frustrated triangular tube, N is 
the number of sites) and the second one enumerates the 
position of the site within the cell, see Fig. [TJ 

The one-particle (one-magnon or one-electron) energy 
spectra for both models (with h — or fi — 0) can easily 
be calculated yielding 

3 

£1,2 («) = ~o^2 - Ji = -£, 



6 3A (k) = -2Jj T Jiy 1 + 3(1 + cosk) (2.3) 

(double-tetrahedra chain) and 
3 

£1,2 («) = ~2^ 2 ~~ 3Jl = ~ £ ' ! 

£3(«) = -3 J\ + 3 Ji cos k (2.4) 
(frustrated triangular tube) for the spin model and 

£ 1,2( K ) = — ^2 = — e i 

£ 3 ,4(k) =*2T^/*2+6*?(1 + cosk) (2.5) 
(double-tetrahedra chain) and 

£1,2 («) = ~te = —e, 
£3(k) = 2^2 + 6h cosk (2.6) 

(frustrated triangular tube) for the electron model. The 
flat (dispersionless) bands £1,2 (ft) = — £ allow to con- 
struct such wave packets of Bloch states which are local- 
ized on the triangles. These localized one-particle states 
read 

\+)rn = ^ (s~, X + UJS m 2 + U} 2 S m ^ 3 ) |FM), 

h)m = ^= + w 2 s- >2 + ws- >3 ) |FM), (2.7) 

where |FM) denotes the ferromagnetic background for 
the spin model and 

\+)m = -7= (c f m l + Ucl l2 + UJ 2 cl n 3 ^ |0), 

\-)m = -J= (4,! + w 2 cL, 2 + "4, 3 ) |0>, (2.8) 
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where |0) denotes the vacuum state for the electron model 
and the spin index a is omitted as irrelevant for the one- 
electron problem. Here u = e 2rrl ^ 3 . 

The two-fold degeneracy of the flat bands corresponds 
to two possible values of the chirality of the triangle. For 
the spin model, after introducing the chirality operator 
for a triangle^ 

4 

Xm = ( s m,l • [s m ,2 x S m ,3]) 

= [\ s m,l s m,2 ~ s m,l s m,2J s m,3 

+ { S m,2 S m,3 — S m,,2 S m,3f S m,l 
+ { S m,3 S m,l ~ S m,3 S m,l) S m,2\ ' (2-9) 

we find Xm|±)m = ±|±)m- We notice that the s z oper- 
ators in Eq. (|2.9[) yield simply 1/2 after acting of \m on 
the states |±) m (|2.7[) . Therefore we may choose a simpler 
form of the chirality operator omitting the operators s z 
and the factor 2 in the last expression in Eq. (|2.9[) . see, 
e.g., Ref. @. Similarly, for the electron models 

if] t 

Xm — ~ \ C mi i c m,2 + c m,lC m 2 

~^ C m,2 C m,3 + c m,2cJ n3 + C^^Cn^i + Cm,Z<^ m ^i (2.10) 

and again Xm|±)m = i|i)m- 111 both spin and electron 
cases the chirality operator \m can be written in the form 

Xm = \+)m(+[,n ~ RmMm, (2.11) 

sec Eqs. ([2/7]) and (|2l)l . 

From the above equations for the spectra it is obvious 
that the two-fold degenerate flat band becomes the lowest 
one if J2 > 2Ji for the spin model or > 2t\ for the 
electron model. In what follows we assume that these 
ratios arc fulfilled. 



III. THE CONTRIBUTION OF INDEPENDENT 
LOCALIZED STATES TO THERMODYNAMIC 
QUANTITIES 

The spin Hamiltonian (|2.1[) commutes with S z , i.e., 
the number of magnons n = N/2 — S z is a good quan- 
tum number. Similarly, the electron Hamiltonian (|2.2[) 
commutes with the operator of the number of electrons. 
Therefore, we may consider the subspaces with different 
numbers of magnons or electrons separately. Moreover, 
we may assume at first h = or /i = and add triv- 
ial contributions of these terms to the partition function 
later. 

We start with the construction of localized many- 
particle eigenstates in the subspaces with n < n max oc TV 
magnons or electrons based on the localized one-particle 
states. These localized many-particle states are obtained 
by occupying the triangular traps with localized parti- 
cles. For the occupation of the traps certain rules have 



to be fulfilled, cf. Ref. for spin systems and Ref. H for 
electron systems. For the spin system on the frustrated- 
tube lattice localized magnons cannot occupy neighbor- 
ing triangular traps whereas for the double-tetrahedra 
spin chain the occupation of neighboring triangular traps 
is allowed. Hence, according to Ref. 0, the frustrated tri- 
angular spin tube belongs to the hard-dimer class and 
the double-tetrahedra spin chain belongs to the hard- 
monomer class. For both electron models localized elec- 
trons may occupy neighboring traps. Moreover, for the 
electron system it is possible that two electrons forming a 
spin-1 triplet (e.g., two spin-up electrons) but having dif- 
ferent chiralities occupy the same triangular trap. Note 
that the different trap occupation rules lead finally to the 
different relations between the maximum number of lo- 
calized magnons (electrons) n max and the number of cells 
Af given below. 

It is helpful to bear in mind a simple picture visualizing 
this construction of the many-particle states^ - — Namely, 
the construction of the many-particle states may be as- 
sociated with a filling of an auxiliary lattice (a simple 
chain of Af sites in all cases considered here) by hard-core 
objects (hard monomers or hard dimers) of two colors 
corresponding to two values of the chirality. Moreover, 
for the electron systems we have to take into account 
in addition the electron spin and the Pauli principle. 
Thus, the maximum filling with magnons is n max = Af 
(double-tetrahedra Heisenberg chain) and n max = Af/2 
(frustrated Heisenberg triangular tube), whereas for the 
Hubbard model the maximum filling with electrons is 
"■max = 2Af for both lattices. 

According to these rules the localized many-particle 
states are product states of localized one-particle states 
with the energy Efm — ne {E~fm is the energy of the fer- 
romagnetic state) for the spin models or with the energy 
— ne for the electron models, where e is given in Eqs. (|2.3[) 

(12. 6p . Importantly, localized electron eigenstates con- 
structed in this way do not feel the Hubbard interaction 
U. Furthermore, the localized many-particle states are 
the only ground states in the corresponding subspaces 
with up to n ma x magnons or electrons if J2 > 2J\ (spin 
models) or i 2 > 2t\ (electron models). We have checked 
this analyzing full diagonalization data for several finite 
spin and electron systems. Obviously, there is a large 
manifold of degenerate localized many-particle ground 
states in an n-particle subspace. We will denote this 
ground-state degeneracy in what follows as gj^{n). 

For the spin models the contribution of these localized 
eigenstates to the partition function is given by 

"■max N 
J-^ X E FM -q-h~ n (h 1 -h) 

Z lm (T,h,N) = 9Mn)e * 

71=0 

JV n max 

= e t 2_^ gjv (n)z n , z = e~^, hi = £,(3.1) 

71 = 

where the quantity g^f{n) represents the degeneracy of 
the ground-state manifold of n magnons in a system with 
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Af traps. For the spin models it is easy to obtain (see 
Ref.i) 



g*(n) = 2 n C% , C% 



Af\ 



n\(Af-n)\ 
(double-tetrahedra chain) and 

g/*(n) = 2 n Z hA (n,M) 



(3.2) 



(3.3) 



(frustrated triangular tube), where Z^^n^Af) is the 
canonical partition function of n hard dimers on a pe- 
riodic chain of Af sites. The factor 2™ in the above ex- 
pressions stems from the extra degeneracy due to the 
chirality degrees of freedom. After substitution of g\r{n) 
from Eq. (pT2]) or Eq. (|3~5|) into Eq. ([3TT]) one obtains the 
free energy 



F lm {T,h,N) E 



FM 



TV 

Af Af W 

(doublc-tctrahcdra chain) or 



h-T]n(l + 2z) (3.4) 



F lm (T, h, N) _ E FM N_ h _ r ln (Af + Af ) 



Af 



Af 2Af 

i /r 



Af 



Al ' 2= 2 ± Vl + 2Z 



(3.5) 



(frustrated triangular tube). At low temperatures and 
for magnetic fields h around the saturation field h\ = e 
the contribution of localized states is dominating. Hence, 
F lm (T,h,N) given in Eqs. ([3~4]) and (|3~5|) yields a good 
description of the low-temperature physics near the sat- 
uration field of the full spin model. 

Analogously, the contribution of the localized eigen- 
states to the grand-canonical partition function of the 
electron models is given by 



n=0 



(3.6) 



n=0 



To avoid the calculation of gjy(n) one may rewrite 
Eq. ()3.6[) as a sum over occupation numbers of each cell 
taking into account (i) that the cells are independent and 
(ii) that each cell may contain 0, 1, or 2 electrons having 
the degeneracy of the ground statesoi(O) = 1, gi(l) = 4, 
or 5i(2) = 3, respectively, see Ref. Thus, we have 

B u (T,n,N) = 

£ ... 9i{nx)...9i{nM)z^ + - +nj " 



E 9i(n)z n 



."=0,1,2 



Af 



= (l + 4z + 3z 2 f (3.7) 



for both lattices, the double-tetrahedra chain and the 
frustrated triangular tube. Eq. (|3.7[) immediately yields 
the required grand-thermodynamical potential 



a c (?>,A0 
Af 



= -Tlnfl +4z + 3z 2 



(3.8) 



Again, at low temperatures and for chemical potentials 
/i around /iq — e the contribution of localized states is 
dominating, and Q\ e (T, /i, N) given in Eq. (|3.8[) yields a 
good description of the low-temperature physics of the 
full electron model. 

We mention that the obtained formulas for the free 
energy and the grand-thermodynamical potential are 
similar (but not identical) to those derived in previ- 
ous papers^ the deviations from the previously derived 
equations are related to the chirality degrees of freedom. 

Let us briefly discuss the low-temperature thermody- 
namics as it follows from Eqs. (|3.4p. (|3.5p . and (|3.8p . The 
main low-temperature features of the spin models for h 
around h\ are as follows: (i) a jump in the ground-state 
magnetization curve at the saturation field hi with a pre- 
ceding wide plateau, (ii) a nonzero residual ground-state 
entropy at the saturation field hi, (iii) a low-temperature 
peak in the specific heat, which moves to T = as h ap- 
proaches hi. Correspondingly, for the electron models we 
have: (i) a zero-temperature jump in the averaged num- 
ber of electrons as a function of the chemical potential 
at /i = fj,Q, (ii) a nonzero residual ground-state entropy 
for (j, = /xo (or as a function of the electron concentration 
c = n/Af for c < 2), (iii) a low-temperature peak in the 
grand-canonical specific heat C(T, /i, A), but a vanishing 
low-temperature canonical specific heat C(T n, N) = 
for n < n max (see also the discussion in Ref. Wh . 

The temperature dependence of the specific heat of 
the spin models is shown in Fig. [2j By comparison of 
the hard-core data obtained from Eqs. (|3.4p . (|3.5p with 
exact-diagonalization data of the full spin model we can 
estimate the range of validity of the hard-core descrip- 
tion. The low-temperature parts of the curves for the 
specific heat obtained from exact diagonalization (sym- 
bols) and from the hard-core description (dashed lines) 
coincide at least up to T — 0.15. For both Heiscnbcrg 
chains the specific heat shows in addition to the typi- 
cal high-temperature maximum around T ~ J2 a low- 
temperature maximum which is well described by the 
hard-core models. This low-temperature maximum can 
be ascribed to an extra low-energy scale set by the lo- 
calized cigenstates. Interestingly, for the frustrated tri- 
angular tube there is even a third maximum which can 
be related to a third energy scale set by another class 
of highly degenerate eigenstates, the so-called interact- 
ing localized-magnon states, which will be discussed in 
the next section. 

There arc no finite-size effects in the hard-monomer 
description (|3.4p . To illustrate the finite-size dependence 
inherent in the hard-dimer description (|3.5[) , we compare 
in panels (c) and (d) of Fig. [5] results for finite Af (dashed 
lines) with data for Af 00 (thin dashed lines). Clearly, 
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finite-size effects arc obvious only at low temperatures 
for h below the saturation field hi. 

The high degeneracy of the localized eigenstates leads 
to a residual ground-state entropy given by S(T = 0, h = 
hi,N)/M = In 3 ps 1.099 (double-tetrahedra spin chain) 
and S(T = 0,h = h u N)/Af = ln2 ps 0.693 (frustrated 
triangular spin tube). Due to the chirality these num- 
bers exceed the corresponding results reported in Ref. U> 
for the standard diamond spin chain S(T = 0, ft. = 
h±,N)/Af = In 2 ps 0.693 and the frustrated two- leg spin 
ladder S(T = 0,h= hi,N)/M = ln[(l + VE)/2] ps 0.481. 

For the Hubbard model we have calculated curves for 
the grand-canonical specific heat similar to those for 
the Heisenberg model, which for the sake of brevity are 
not shown here. They also exhibit an additional low- 
temperature maximum for fj, < and M ^ Mo which is 
well described by the localized eigenstates, see also Ref.H. 
The residual ground-state entropy lim.T->o S(T, n, N) = 
In gj\f(n) as a function of the electron concentration c = 
n/Af for the Hubbard model is shown in Fig. [3] Accord- 
ing to Eq. (|3.6[) for finite systems one has n\g_\f(n) = 
d n E ie (z,N)/dz n \ z = . Using Eq. ([377]) one obtains flw(n) 
for (small) finite Af. For Af = 4 the analytical predic- 
tions which follow from Eqs. (|3.7[) . (|3.8|) coincide with 
exact-diagonalization data for the full Hubbard model, 
e.g., g 4 (n) = 1, 16, 108, 400, 886, 1200, 972, 432, 81 for n = 
0, 1, 2, 3, 4, 5, 6, 7, 8. Using Eq. (|3.8[) by means of standard 
relations of statistical mechanics we find for the residual 
ground-state entropy in the thermodynamic limit S(T = 
0,n,N)/Af = ln(l + 4z + 3z 2 ) - (Azlnz + 6z 2 lnz)/(l + 
4z + 3z 2 ) with z = [2(1 - c) - Vc 2 - 2c + 4]/[3(c - 2)], 
c = n/Af, see the dotted curve in Fig. [3] This quantity 
reaches its maximum, 3 In 2 ps 2.079, at c = n/Af = 5/4. 
Moreover, it equals In 3 ps 1.099 at c = n/Af = 2. Ob- 
viously, the chirality leads to an increase of the resid- 
ual ground-state entropy. While for the Hubbard model 
on the diamond chain and the frustrated two-leg ladder 
S(T = 0,n = fj, , N)/N = In 3 ps 1.099 (see Ref. lb), for 
the considered lattices S(T = 0, u = u , N) /Af = 3 In 2 ps 
2.079. 



IV. BEYOND INDEPENDENT LOCALIZED 
MAGNONS 

For the frustrated-tube spin model we can extend the 
hard-dimer description taking into account additional 
low-energy states following the lines described in Ref. [22|. 
If one allows the occupation of neighboring traps by lo- 
calized magnons (i.e., relaxing the hard-dimer rule) one 
has also an eigenstate of the spin Hamiltonian, how- 
ever with a higher energy. More precisely, if two lo- 
calized magnons become neighbors the energy increases 
by J\. Importantly, these localized-magnon states are 
also highly degenerate and they are the lowest excita- 
tions above the independent localized-magnon ground 
states for S z = N/2, N/2 - Af/2 if J 2 > J| (strong- 
coupling regime). Based on finite-size calculations for 
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FIG. 2: (Color online) Specific heat C{T,h,N)/Af vs tem- 
perature T for h around hi for [(a) and (b)] the N — 20 
double-tetrahedra spin chain and [(c) and (d)] the N — 18 
frustrated triangular spin tube. Ji = 1, J% = 5; symbols 
correspond to exact-diagonalization data, dashed lines corre- 
spond to independent localized-magnon predictions derived 
from Eqs. (|3.4|l and (|3 . 5 1) (abbreviation lm), dotted lines cor- 
respond to interacting localized-magnon predictions derived 
form Eq. (|4.2p (abbreviation LM, for the frustrated triangular 
tube only). The thin dashed and thin dotted lines in panels 
(c) and (d) correspond to localized-magnon predictions in the 
limit Af — > co. Note that some curves in panels (c) and (d) 
practically coincide. Namely in panel (c), thin dashed and 
thin dotted lines are indistinguishable at low temperatures, 
whereas at higher temperatures dashed and thin dashed (dot- 
ted and thin dotted) lines coincide. In panel (d) dashed and 
thin dashed (dotted and thin dotted) lines cannot be distin- 
guished. 



N = 18,..., 72 we estimate Jf/Ji ps 2.68 > 2. These 
additional eigenstates can be described as interacting 
localized-magnon states, where the repulsive interaction 
V = Ji is responsible for the energy increase with respect 
to the independent localized-magnon states. Taking into 
account this finite repulsion V in the partition function 
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FIG. 3: (Color online) Residual ground-state entropy S(T = 
Q,n,N)/N vs electron concentration c = n/Af for the Hub- 
bard model with t2 > 2ti on the double-tetrahedra chain and 
the frustrated triangular tube (the data for both systems are 
identical). Note that the curves becomes indistinguishable in 
the scale of the figure as M increases. 
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FIG. 4: (Color online) Specific heat C(T, h, N)/N vs temper- 
ature T for finite frustrated triangular spin tubes of N = 18 
sites (M = 6 cells) with Ji = 1, J 2 = 5 {hi = 10.5) for 
various magnetic fields down to values below /12 = 8.5 (lines 
- lattice-gas model with finite repulsion, symbols - exact- 
diagonalization data for the full spin model). 



of the lattice-gas model with nearest- neighbor interaction 
we have then [instead of Eq. (|3.1[) ] 

E - N h 

Z LM (T,h,N) = e-^T^ X) E gi ^ 

Tli— 0,1 n M— Oil 

/ \ 77,-1- -4-n.r V(n 1 Ti 2 + Ti 2 n 3 + ... + Ti A rTi 1 ) 

x • • .gi{n^)z ni+ - +nu e * (4.1) 

with 31(0) = 1, 51(1) = 2 and z is defined in Eq. ([3d]). 
In Eq. (|4.1[) we have used a representation in terms of 
the cell occupation numbers rii, cf. Eq. (|3.7|) . Obviously, 
Eq. (|3.1|) is obtained from Eq. (|4.ip for y — )• 00. Evaluat- 
ing the sums in Eq. (|4. 1[) by means of the transfer-matrix 
method we arrive at the following result for the free en- 
ergy 

F LM (T,h,N) _E FM N ln(Af + Af ) 
Af N Uf_ 7V" 

Ai,2 = ^ + ze~T ± </i + 2z — ze'T 1 + z 2 e~ ? r .{4.2) 

Note that the lattice-gas model with finite repulsion 
fi~ll , (|4T2"]) takes into account 3^ ps 1.442^ states^ 
of the 2 W eigcnstates of the initial quantum spin model 
(12. ip . whereas hard-dimer model (|3.1[) . (|3.3p . (|3.5|) has 
only 2^ ps 1.260^ states^ On the other hand, the 
hard-monomer model ([3TT]) . ([3"^]) . (JSH) has 3^ ps 1.316^ 
states^ 

The specific heat derived from Eq. (|4.2[) is plotted in 
Fig. [3J panels (c) and (d), dotted and thin dotted lines. 
Indeed, the inclusion of the interacting localized-magnon 
states leads to a significant improvement of the lattice-gas 
description. The lattice-gas model with finite repulsion 
(|4.1[) , (|4.2p covers the thermodynamics of the full spin 
model at least up to T = 0.9 for h ~ hi including the 
two maxima below the typical high-temperature maxi- 
mum around T ~ Ji. Again finite-size effects are more 
important at low temperatures for h below the satura- 
tion field hi, see thin dotted lines in panels (c) and (d) 

of Fig. m 



The interacting localized-magnon states being excita- 
tions for S z = N/2, . . . , N/2 — Af/2 can become ground 
states for smaller values of S z . The ground-state mag- 
netization curve (M) = 2S Z /N vs h for the frustrated 
triangular spin tube presented in Ref. [ID exhibits two 
plateaus at (M) = 2/3 for h 2 < h < hi and at (M) = 1/3 
for h^ < h < hi and two jumps at h = hi and h = hi, 
where hi = 3Ji + 3J2/2, h 2 = Ji + 2>J%/2, and h 3 = 2Ji. 
In the language of localized magnons the 2/3-plateau 
corresponds to the maximum filling with n = n max = 
TV/ 2 independent localized magnons (i.e., every second 
trap is occupied and S z = N/2 - Af/2 = N/3). If 
n max <n< 2n max (i.e., N/2-N < S z < N/2 - Af/2), 
the ground states in the strong-coupling regime are ob- 
tained by filling the remaining empty cells (i.e., the hard- 
dimer rule is relaxed) thus having interacting localized- 
magnon states as ground states. Then the very broad 
1 / 3-plateau corresponds to the complete filling of all cells 
with n = 2n max = TV localized magnons. Hence, with 
the improved effective theory given in Eq. (|4.2p we can 
provide an accurate description of the low-temperature 
physics of the frustrated triangular spin tube in the 
strong-coupling regime not only near the saturation field 
hi up to quite large temperatures as shown in Fig. [2 but 
also for much lower magnetic fields in the entire region 
of the 2/3-plateau and even for fields within the 1/3- 
plateau being not too far from h 2 , see Fig. 5] It might be 
interesting to recall that the lattice-gas model with finite 
repulsion provides similar description of the frustrated 
two-leg spin ladder around both characteristic fields hi 
and h 2 *£ This is not the case for the frustrated triangu- 
lar spin tube because of the chirality degrees of freedom, 
compare the results for h = 10 and h = 9 in Fig. 0J 



V. LIFTING THE DEGENERACY DUE TO 
CHIRALITY DEGREES OF FREEDOM 

As discussed above the chirality degrees of freedom 
lead to an extra degeneracy of the independent localized- 
magnon or localized-electron ground states in the sub- 
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space with n < n max magnons or electrons. Moreover, 
the interacting localized-magnon states discussed for the 
frustrated triangular spin tube also carry this extra de- 
gree of freedom. This degeneracy of the eigenstates ow- 
ing to the chirality may be lifted by a small symmetry- 
breaking perturbation. As a rule, perturbations of ideal 
model Hamiltonians may lead to a more realistic descrip- 
tion of real systems. We consider here separately for 
the spin systems the case of a Zeeman-like perturbation 
(which corresponds to a Dzyaloshinskii-Moriya interac- 
tion between neighboring spins in a triangular trap), see 
Sec. IV Al and the case of an XX-like perturbation acting 
on (pseudo)spin variables representing chiralities (which 
corresponds to a four-site interaction between spin pairs 
in neighboring traps), see Sec. IV CI For the Hubbard 
model appropriate perturbations correspond to a mag- 
netic field perpendicular to the triangular traps and a 
four-site electron-electron interaction, see Sec. IV Bl and 
Sec.lVCl 



A. Dzyaloshinskii-Moriya interaction 

We consider the spin system in the subspaces with 
n < n max magnons. The ground state has the energy 
Efm — ne (h = 0) and the degeneracies c/aK 71 ) are given 
in Eqs. (pH|) or (f3~3")l . The factor 2™ present in these 
formulas for gj^(n) is due to the chirality of localized 
magnons. We introduce the (pseudo)spin-l/2 operators 

< n = 2 (l+)m(+|m - RmHm) = \x m , (5.1) 

where Xm = (t/VsKs^s^ - s~ nA s^ 2 + s+ >2 s^ 3 - 

s m,2 s m,3 + S m,3 s m,i ~ S m,3 s m,i) is the chirality operator, 
see Eq. (|2.9p and the discussion below this equation. We 
now add to the spin Hamiltonian (|2.1[) a small perturba- 
tion 



2V3 



/ ^ { S m,l S m,2 s m,l 5 m,2 



S m,2 S m,3 S m,2 S m,3 + S m,3 S m,l S m.3 S m,l) 
) V (i x s y - s y s x 



-a x q y _ q y q x 4_ o 1 <?y _ <?y \ 

s m.2 s m,3 b m,2 b m,3 ' b m,3 b m,l 'm.sVli ' 



(5.2) 



where the last expression in Eq. (|5.2|) corresponds to 
a Dzyaloshinskii-Moriya interaction D = De z , D = 
eW/\/3 between neighboring spins within each triangu- 
lar trap. Note that the perturbation commutes with 
S z . According to the discussion of the chirality operator 
(|2.9p in Sec. [TT] it is obvious that the localized states are 
also eigenstates of the perturbation Hamiltonian (|5 . 2|) . 
The set of 2™ degenerate ground states in the unper- 
turbed system belonging to one particular spatial config- 
uration of n magnons placed in a certain (allowed) set of 



traps now splits into n + 1 subsets of levels. The subsets 
are characterized by the magnon numbers n+ and ro_ 
(n + + Ti- = n), belonging to the chirality indeces + and 
— , respectively. There are n!/(n + !?i_!) degenerate states 
in the subset with energy Efm— ne+(n+— n^)e^ /2. The 
effective Hamiltonian acting in the subspace of the gjy(n) 
former n-magnon ground states of the unperturbed sys- 
tem reads 



n s + = E] 



FM 



1} E- 



(5.3) 



where the sum runs over the n occupied traps only. 

We consider now the partition function of the spin 
model with the Hamiltonian H s + H^, Eqs. (|2.1[) 
and (|5.2p . at low temperatures and h close to hi- 
The dominant contribution to the partition function 
of the spin system comes from the low-energy de- 
grees of freedom, which are governed by the Hamil- 
tonian (|5.3|) . Therefore the partition function is 
given by Eq. (|3.1|) replacing gx{n) by gtf{ri) = 
{2cosh[eW/(2T)]} n Q (double-tetrahedra chain) and 
gM{n) = {2cosh[e( 1 V(2T)]}"Z hd (n,AA) (frustrated tri- 
angular tube), cf. Eqs. (|3.2p and (|3.3p . As a result, the 
free energy reads 



F lm (T,h,N) E FM 



N 



A'' 



M 2Af 



h — Tin 1 



e (i) 

2 cosh— z$A) 



(double-tetrahedra chain) and 



F lm (T,h,N) 



E, 



FM 



N h r ln(Af + A^) 



A' 



A" 2Af 



Ai,: 



1 



± 



1 n , 

_ + 2 cosh— z 



(5.5) 



(frustrated triangular tube), cf. Eqs. (|3.4p and (|3.5p . 
Similar to Sec. IIVI wc can take into account inter- 
acting localized-magnon states also for the perturbed 
frustrated triangular tube described by the Hamiltonian 
H s + H^K The improved free energy -Flm(7\ h, N) is 
then given by Eq. (|4.2p . where 2z has to be substituted 
by 2cosh[eW/(2T)]z. 

In Fig. [5] we compare exact-diagonalization data 
for perturbed spin systems with predictions based on 
Eq. (|5.4p and improved Eq. (|5.5p . Small perturbations 
lead to splitting of the ground-state levels and therefore 
to arising of one more low-energy scale. As a result, low- 
temperature features close to hi are more subtle. Thus 
temperature profiles of the specific heat for the spin sys- 
tems show more tiny features which can be seen in Fig. [5] 
For a special set of parameters the temperature depen- 
dence C vs T may exhibit even three (four) maxima for 
the double-tetrahedra chain (frustrated triangular tube). 
The low-temperature maxima arc excellently described 
within the effective low-energy theory, see Fig. [5] 
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FIG. 5: (Color online) Specific heat C(T, h, N) /Af vs temper- 
ature T for ft = 0.98/ii for the spin model (|2.1j) . (|5.2|) on (a) 
the double-tetrahedra chain and (b) the frustrated triangular 
tube. Exact-diagonalization data (symbols) versus analytical 
predictions (lines) according to Eq. (|5.4[) (double-tetrahedra 
chain) and improved Eq. (|5.5[) (frustrated triangular tube) 
for finite systems of N = 4 cells with Ji = 1, J% = 5, and 
D = 0, 0.025, 0.05, 0.1. 



B. Electrons in a magnetic field 

In analogy to the above discussion for perturbed spin 
models, we consider the case of n < n max Hubbard elec- 
trons (|2.2p with a perturbation 

<T=t4 

^Q 1 ] = 'jX! ( C L,1,,t C ™,2,o- + C m ,l,ir4,2, l T 
m 

+ c m,3,o- c rra>l)CT + Cm.S^m.ljO- ) i (5-6) 

where is a pure imaginary component of the hopping 
integral between neighboring sites along the triangular 
traps. For the perturbed Hamiltonian the number of 
electrons remains a conserved quantity and the localized 
states are its eigenstates, cf. the chirality operator (|2.10p 
in SecHU The perturbation considered here corresponds 
to a magnetic field perpendicular to the triangular trap. 
Then the hoping integral ty acquires the Peicrls phase 
factor e ( 2 W*o) /, rfr-A^ w k ere <j) _ ^ c / e j s g ux q Uan _ 

turn, A is the vector potential of the external magnetic 
field, see Refs. I25H2E The effective Hamiltonian acting 
in the subspace of the gj^{n) former n-electron ground 
states of the unperturbed system reads 

We + W^-ne + eW^i (5.7) 

m 



where e*- 1 ) = —2\/Zg and the sum runs over the n occu- 
pied traps only. 

We consider now the grand-canonical partition func- 
tion of the electron models with the Hamiltonian H c + 
H^, Eqs. (|2.2|) and (|5.6[) . at low temperatures and /x 
close to no = e. The dominant contribution to the grand- 
canonical partition function comes from the low-energy 
states, which are governed by the Hamiltonian (|5.7[> . Re- 
peating the arguments which lead to Eq. (|3.8|) we arrive 
now at 

4M-Th( 1+ «^ + tf), M 

where z = eS £ ~^/ T , cf. Eqs. (j3H) - (03]). 

We focus on the low-temperature behavior of 
the entropy of the perturbed electron model. 
We can easily find, using a simple counting of 
states, the residual ground-state entropy, namely, 
S(T = 0,™, AT) = ln(2"C^) for n = 0, . . . ,W and 
S(T = 0,n,A0 = ln(3 n - Ar 2 aV - n C^ _ ' V ) for n = 
Af, . . . , 2Af. In the thermodynamic limit Af — > 00 this 
gives: S(T = 0,n,Af)/Af = cln2-clnc-(l-c)ln(l-c) 
for < c < 1 and S(T = 0,n,Af)/Af = 
(c-l)ln3 + (2-c) ln2-(c-l)ln(c-l)-(2-c) ln(2-c) 
for 1 < c < 2, where c = n/Af. For the spe- 
cial electron concentrations c = 1 and c = 2 
one has S(T = 0)/Af = In 2 w 0.693 and 
S(T = 0)/Af = In 3 w 1.099, respectively. For finite tem- 
peratures we find from Eq. (f578|) S(T, n,Af)/Af = ln{l + 
4cosh[e( 1 ) /(2T)]z + 3z 2 } - {2(e« /T) smh[e« /(2T)]z + 
4 coshfeW /(2T)]z In z + 6z 2 In z} /{1+4 cosh[eW /{2T)]z + 
3z 2 } with z = {2(1 - c) cosh[eW /(2T)] - 

^4(c - l) 2 cosh 2 [e( 1 )/(2T)] - 3c(c - 2)}/[3(c - 2)]. 

In Fig. [6] we show the entropy S(T,n,N) of the per- 
turbed electron systems obtained by the above given for- 
mulas versus the electron concentration c = n/Af. To es- 
timate the region of validity of these results we compared 
first exact-diagonalization data for the grand-canonical 
specific heat C(T, a, N) for finite perturbed electron sys- 
tems (e.g., for Af = 2, ti = 1, <2 =3 and ti = 5, 
g = 0.1, fi = 0.98/io, Mo, 1-02/Zo) with analytical results 
for C(T, fj,, N) based on Eq. (f5T8]> . (For the sake of brevity 
we do not show these results explicitely.) We find an ex- 
cellent agreement between these results at least up to 
T = 0.2 for ti = 3 (and even for higher temperatures 
for larger values of t-^). For nonzero but low tempera- 
tures (e.g., T = 0.15) S(T, n,Af)/Af behaves as in Fig. El 
see panel (c) of Fig. [5] However, at lower temperatures 
(e.g., T = 0.075 and T = 0.015) the smallest energy 
scale comes into play, and the entropy changes remain- 
ing nonzero in the ground state, see panels (b) and (a) 
of Fig. [5] In spite of the partial degeneracy lifting due to 
the perturbation (|5.6p . the ground states remain hugely 
degenerate and exhibit a nonzero residual ground-state 
entropy for electron concentrations < c = n/Af < 2, 
see panel (a) of Fig. © (T = 0.015). 
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nian (12.11) reads 



(b) 2 



1? 

B 
to 



0.5 



7=0.075 






1 

r — 0' ^ 








5^32 




5^16 ~A— 




3V^=8 — * — 



0.5 



(c) 2 
1.5 
1 



1 

nisi 



1.5 



0.5 




0.5 



1 



1.5 



r=o.i5 




'X 










5^32 




9\£=16 ~A— 




5\£=8 — * — 



FIG. 6: (Color online) Entropy S(T,n, Af) /N vs electron con- 
centration c = n/Af at low temperatures (a) T = 0.015, (b) 
T = 0.075, and (c) T = 0.15 for the Hubbard model (|23j) 
with t2 > 2ti completed by small perturbation (|5.6[) with 
j = 0.1 on the double-tetrahedra chain and the frustrated tri- 
angular tube. Analytical predictions are obtained according 
to Eq. (|5.8|) for systems of Af = 8, 16, 32, and Af — > oo cells. 



e (2) 

ff s (2) = -9" 51 [( S ™,1 + WS ™,2 + ^ 2 s m , 3 ) 

rn 

X ( S m,l + WS m,2 + ^m.a) 
X ( s m+l,l + w2s m+l,2 + U}S m+l,a) 

x + w 2 s+ +1)2 + ws+ +li3 ) + H.C.] , (5.10) 

where the sum runs over all trapping cells and w = e 2 ™/ 3 . 
Obviously, that corresponds to certain four-site interac- 
tions with the interaction constant 6^/9. The perturba- 

(2) 

tion Hamiltonian Hs commutes with S z . Moreover, af- 
ter acting on the localized-magnon crystal state the per- 

turbation Hamiltonian Hs changes the chirality indeces 
only. The effective Hamiltonian acting in the subspace 
of the localized-magnon crystal states of the unperturbed 
system now reads 



H s + Hi 2) = E- 



FM 



A/" 

,(2) \p I + - 



H.C.I5.11) 



Due to the perturbation H s 2 ^ the 2-^-fold degenerate 
ground state of the double-tetrahedra spin chain with 
n = n max = Af magnons splits into Af + 1 groups 
of sublevels. Moreover, the effective (pseudo)spin-l/2 
XX chain (|5.1ip is the exactly solvable model 2 ^ and 
therefore we immediately obtain for the partition func- 
tion Z(T,n,N), n = N/2 - S z of the double-tetrahedra 
spin chain with a Hamiltonian given by the sum of the 
terms in Eqs. (|2.1[) and (|5.10[) the following dominant 
contribution at low temperatures for the magnetization 
S z = N/2 - rw = N/4 



Zi m (T,n max ,7Y) = e~ 



~[ 2 cosh 



v- 2 ' cos k 
f 



(5.12) 



C. Interacting pseudospins 

Now we consider the double-tetrahedra spin chain at 
the (M) = 1/2 plateau (i.e., there are n = n max = 
Af magnons in the spin system, that is the so-called 
localized-magnon crystal state). We add a perturbation 
which can be understood as an XX interaction of chi- 
rality pseudospins on neighboring trapping cells. The 
(pseudo)spin raising and lowering operators are defined 
as 



(5.9) 



They can be expressed by bilinear forms in the spin op- 
erators s~ and s + attached to the mth cell, see Eq. (|2.7p . 



Clearly, r+|+) 

7 



T H 

"5 1 n 



[|+)m = i->, 



K = 2irl/N, I = -N/2, -N/2 + 1, . . . ,N/2 - 1 (we as- 
sume without loss of generality that N is even). Low- 
temperature thermodynamic quantities in the thermo- 
dynamic limit N — > 00 follow from the free energy 



± lm 

(T AO 
N 



E 



FM 



N 
cosk 



■ — / dn In I 2 cosh ■ 



0. The perturbation added to the Hamilto- 



(5.13) 



Note that the perturbation (|5.2j) can also be included. 
Then the formula (|5.13|) has to be slightly modified, 
namely, 2e' 2 ' cos k in Eq. (|5 . 13[) has to be replaced by 
e' 1 ) + 2e( 2 ' cos k. Thus, including both perturbations 
(|5.2[) and (|5.10|) we arrive at an effective (pseudo)spin- 
1/2 XX chain in a transverse field which governs the 
low-temperature physics of the double-tetrahedra chain 
in the subspace with n = ro max = /V magnons (i.e., at 
the magnetization S z = N/2 — n max = N/4). It might 
be interesting to mention here that the spin- 1/2 XX 
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chain in a transverse field also emerges as an effective 
low-energy model for the diamond spin chain at high 
magnetic fields if the conditions ensuring the presence of 
localized magnons become slightly violated^ For such a 
generalized diamond chain the spin-1/2 transverse XX 
chain describes a weak spreading of the independent lo- 
calized magnons over the whole chain. Another related 
model, the so-called spin-chirality model, is used for ef- 
fective description of three-leg spin tubes within the per- 
turbation theory approach from the strong rung-coupling 
limit, see Ref. [l7| and references therein. We stress here 
that in the case at hand the spin-1/2 transverse XX 
chain describes propagating of chirality over the whole 
chain in the localized-magnon crystal state. 

We consider next the frustrated triangular tube again 
in the subspace with n = n max = Af/2 magnons (i.e., 
at the (M) = 2/3 plateau). The ground state, besides 
the degeneracy owing to chirality, is two-fold degener- 
ate, i.e., the independent localized magnons may oc- 
cupy either even-site or odd-site sublattice only. The 
perturbation to lift the degeneracy of the localized- 
magnon crystal state with respect to chirality that has 
to be added to the Hamiltonian (|2.ip corresponds to 

Ui 2) = e (2) Em=i ( r m T m+ 2 + H - c -)> i- e -; & represents 
now an XX interaction between next-nearest-ncighbor 
(pseudo)spins. In the initial spin model it is a certain 
four-site interaction which contains the sites of next- 
nearest-neighbor cells. 

Finally we discuss briefly corresponding perturba- 
tions for the electron systems in the sector of n = Af 
electrons^ To have exactly one electron per cell we in- 
troduce an extra repulsion between electrons in neighbor- 
ing cells (i.e., we consider an extended Hubbard model). 
Since the chirality and spin of electron states in each cell 
are not fixed, we have a 4^-fold degenerate ground state. 
A perturbation that lifts the degeneracy owing to chiral- 
ity (independently of the spin) again corresponds to an 

(2) 

XX interaction between (pscudo)spins given by T-L c = 

E CT =t,^S- W$ - £ (2) E™=i (r+r" +1 + H.c), where 
t± are defined by Eqs. ([5T5]l and (|2~g|) . Note that there 
are no spin indeces in the r.h.s. of the formula for 

(2) 

T-Lq a , i.e., each state of the perturbed Hamiltonian is 
2-^-fold degenerate owing to the electron spin. Thus, if 

(2) 

the perturbation interaction He (which contains, gen- 
erally speaking, four-site terms in the electron-electron 
interaction between the neighboring cells) is switched on, 
the thermodynamic properties of the extended Hubbard 
model with n = N electrons on both considered lat- 
tices are related to those of the (pseudo)spin-l/2 XX 
chain, see the corresponding results for the spin model, 
Eqs. (|5T2|) and (pTTB")) . 

VI. CONCLUSIONS 

To summarize, we have considered the low- 
temperature properties of the spin-1/2 antiferromagnetic 



Heiscnbcrg model and the repulsive Hubbard model 
on two ID lattices containing equilateral triangles. 
The lattices under consideration have a dispersionless 
lowest-energy band for the one-particle problem, and 
the corresponding localized one-particle states can be 
trapped on the triangles. Due to the triangular geometry 
of the trapping cells the localized one-particle states 
are characterized by two possible values of the chirality. 
Using the localized nature of the one-particle states we 
can construct corresponding many-particle low-energy 
states. Moreover, we estimate their contribution to ther- 
modynamics exploiting classical lattice-gas description 
of the low-energy degrees of freedom of the quantum 
models. The lattice-gas description yields explicit 
analytical formulas for thermodynamic quantities at 
low temperatures in a certain region of the magnetic 
field (chemical potential) for the spin (electron) model. 
We investigate the effects of the localized states on the 
low-temperature thermodynamics. In detail we discuss 
the specific heat C(T,h,N) for the spin systems and 
the entropy S(T,n,N) = J Q T dT'C(T',n,N)/T' for the 
electron systems. Both quantities exhibit fingerprints 
of highly-degenerate localized states, namely, additional 
low-temperature peaks of C(T, h, N)/N or a finite resid- 
ual ground-state entropy S(T = 0,n,N)/N. Since the 
considered systems show a significant zero-temperature 
entropy, they may exhibit an enhanced magnetocaloric 
effect M 

The degeneracy related to the chirality degrees of free- 
dom may be lifted by small symmetry-breaking inter- 
actions. For the perturbed system we provide an ef- 
fective description of low-energy degrees of freedom of 
the considered spin and electron models in terms of a 
(pseudo)spin-l/2 XX chain in a transverse field. It 
might be interesting to note that in contrast to usual 
cases, where the spins are related to the spin degree of 
freedom of electrons, the (pseudo)spins emerging in our 
case are related to the charge degree of freedom of elec- 
trons and they simply stand for a (pseudo)spin repre- 
sentation of the chirality. Moreover, the chirality inher- 
ent in the considered spin models on geometrically frus- 
trated lattices may also give rise to a chain of quantum 
(pseudo)spins 1/2. 

It is worthy noting that quantum spin chains are of- 
ten used in quantum information theory both for illus- 
tration of basic concepts and as candidates for physical 
implementation^ From such a perspective, manipula- 
tion with chirality^ realized in (pseudo)spin chains may 
be an interesting subject for further studies. Finally, al- 
though the main advantage of the considered strongly 
correlated lattice models is the possibility to elaborate a 
theoretical description of thermodynamics which works 
perfectly well at low temperatures for high magnetic 
fields or low concentrations of electrons, we may mention 
here some experimental solid-state realizations of similar 
systems, see Refs. H6ll33ll34l . 
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